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We investigate the interaction of two Mn ions in the dilute magnetic semiconductor Gai-ajMn^As 
using the variational envelope wave function approach within the framework of six band model of 
the valence band. We find that the effective interaction between the Mn core spins at a typical 
separation d is strongly anisotropic for active Mn concentrations less than Xactivo = 1.3 %, but 
it is almost isotropic for shorter distances (d < 13 A). As a result, in unannealed and strongly 
compensated samples strong frustration effects must be present. We also verify that an effective 
Hamiltonian description can be used in the dilute limit, x act i vc < 1.3 %, and extract the parameters 
of this effective Hamiltonian. 



I. INTRODUCTION 



Emergence of semiconductor with ferromagnetic properties^ in (III-Mn)V materials leads to possibility of inte- 
gration of processing and magnetic storage in a single device. Gai_ x Mn x As is promising material for such devices 
although today Curie temperature up to 140 K is not high enougbA GaMnAs is an example of (III-Mn)V family of 
materials which evolved from II- VI based diluted magnetic semiconductors^. These materials were extensively stud- 
ied in the 80's. Especially from our point of interest, spin orbit interaction was studied in Ref. (see also references 
therein) . Marriage of ferromagnetic and electronic degree of freedom is important today because manipulation of the 
spin in solid state physics starts to be applied at room temperature (recently room temperature fcrromagnctism has 
been observed in GaMnN exhibiting Curie temperature up to 370K reported in Ref.@) with the perspectives to build 
spin diode, spin transistor and eventually quantum computers. 

Ferromagnetic semiconductor like Gai-^Mn^As is characterized by presence of localized magnetic moments and 
simultaneously by presence of mobile holes which mediate ferromagnetic indirect exchange interaction between Mn 
moments 7 . It was shown that substitutional Mn atom (Mn substitute Ga) forms localized magnetic moment of 5/2 
which come from its five 3d-electrons. Additionally, divalent Mn atom substituting trivalent Ga atom introduces one 
hole and forms effective mass acceptor. This idealistic scenarios is not realized in real live and not all intentionally 
incorporated Mn atoms come into Ga sites&. Gai-^Mn^As is highly compensated presumably due to interstitials^ and 
antisite defects (As substituted for Ga) . It means that one to one correspondence between Mn moments and number 
of holes is broken; For typical Manganese concentration x=5% number of holes p=0.3 per Mn atom is expected^. 

There are two natural limits for which simulations are performed: heavilyiiii2*i£ and weakly doped regimesi^*i^. 
In the case of heavy doping one considers valence-band holes moving in disordered potential of defects. In this 
limit Zarand and Janko^Si suggested that frustration effects may be important for the magnetic properties of the 
Gai-xMn^As. Using spherical approximation of the valence band and performing mostly analytical calculations they 
showed that spin anisotropy can suppress magnetization. Calculations made by Brey and Gomez-Santos 11 using full 
6-band valence-band Hamiltonianii show that spin anisotropy is smaller that predicted in Ref. liH In the second limit 
of light doping the local impurity levels overlap only weakly forming impurity bands in the gap. In this case tight 
binding model is usedi£*i£ or local density approximation approach is applied 7 . This limit is also interesting from 
metal-insulator transition point of viewi2ii2i2S. 

There is some controversy in the literature. For example the calculations of Ref. ITI1IT2TIT9I and0 are only valid in 
the high concentration limit of the holes. Furthermore, most previous calculations ignored the very large Coulomb 
potential of the negatively charged Mn impurity, which actually provides the largest energy scale in the problem—' 22 . 
This Coulomb potential can be treated non-perturbatively in the very dilute limit, where an impurity band picture 
appliesi^iSiSSiSiSi. However, the spin-orbit coupling being large, one has to incorporate this also in realistic calcula- 
tions. While the spin-orbit coupling has been completely neglected in Ref. [IfJ Ref. |23 took into account the effects 
of spin-orbit coupling only within the framework of the so-called spherical approximation^, where the spin-orbit 
splitting between the spin j = 3/2 and spin j = 1/2 valence bands is taken to infinity, and therefore anisotropy effects 
are overestimated. 

In this paper we shall investigate the interaction between two Mn ions, using the full six-band Hamiltonian 17 ^ 6 ^ 27 , 
which is thought to give a good approximation for the valence band excitations over a wide energy range. To 
determine the effective interaction, we shall compute the spectrum of 'molecular orbitals' on the two Mn ions through 
the application of variational methods. These molecular orbitals of the MnMn dimers are close analogues of the 
molecular orbitals of the molecule^. As we shall see, the energy of these molecular orbitals depends on the 
separation of the two Mn ions and the direction of their spin, which we treat as classical variables. Knowledge of 
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the spin-dependence of this molecular spectrum allows us to compute the effective interaction between the Mn spins 
provided that there is a single bound hole on the Mn? 'molecule'. 

As we shall also see, both the approaches of Ref. [island Ref. 113 fail: While at large separations the structure of the 
molecular orbitals and the effective interaction is roughly captured by the spherical approximation of Ref. ^3 (giving 
very large anisotropies), for small Mn-Mn separations spin-orbit coupling effects are not substantial, and the Mn-Mn 
interaction is essentially isotropic, though the degeneracy of the bound hole states is different from the one obtained 
through the naive spin 1/2 approach of Ref. 0. We shall also check if it is possible to quantitatively describe the 
spectrum of the Mn2 molecule by an effective interaction, as proposed in Ref. 23. A detailed analysis of the spectrum 
reveals that such an effective model description only works for large Mn-Mn separations, d > 13 A, corresponding to 
an active Mn concentration, £ ac ti V e < 1-3 %■ 

The paper is structured as follows. We will first consider Mn-dimcr 18 problem using envelope wave function 
approach. This approach is justified by the fact that when a small portion of Gallium atoms is substituted randomly 
by manganese atoms we expect that the band model of GaAs will be applicable to Gai-^Mn^As, too. After calculating 
energy states of the dimer we will map them into the effective spin Hamiltonia n 1 ^ 2 ? and discuss results of this mapping. 
We also will discuss spin anisotropy effects. 



II. VARIATIONAL CALCULATIONS 



To investigate the electronic structure of an Mn2 'molecule', we write the Hamiltonian describing the interaction of 
a hole with the dimer as 
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Here the Kohn-Luttinger Hamiltonian Hkl describes the kinetic energy of the holes in the valence band within the 
envelope function approach 2 ^*^, and its detailed form is given in Appendix ^ The Coulomb interaction with the 
negatively charged Mn cores is accounted for by the terms e 2 /eri, with e = 10.67 denoting the dielectric constant 
of GaAs host semiconductor, and r, = |r — Rj,\ being the distance between the two Manganese ions at positions Ri 
(i = 1,2) and the hole at position r. The so-called central cell corrections 2 ^, V c (ri), are used to take into account 
the atomic potential in the vicinity of the Mn core ions. For this central cell correction we used the simple Gaussian 
form, V c (r) = V c e~ r / r ". Finally, the last term in Eq.QJ takes into account the exchange interaction with the two 
Mn core spins, with Si and 5*2 being spin 5/2 operators corresponding to half-filled <i-lcvels. These core spins shall 
be replaced by classical spins in our calculations, and we shall only treat the hole-spin quantum-mechanically. 
To regularize the the exchange interaction in Eq. Q we shall use a slightly smeared delta function, 

0) 

where r p d is a cut-off in the range of the inter-atomic distance. While a microscopic derivation results in a slightly 
more complicated form for exchange interaction^ 2 ^, Eq. © can be also used as long as the spatial scale of the 
exchange interaction, r p d, is small enoughiii^. 

To determine the spectrum of Eq. JQ| we used a variational approach, where we expressed the wave function of the 
holes in the following form, 
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Max) = sj^z-^ 2 ' 2 *^)- 



(4) 



Here the H^s denote Hermit polynomials, and the label /i refers to the six spinor components (see Appendix [SJ. The 
minimization was performed over the parameter a as well as over the linear coefficients, C^ k , k = O...N max — 1). 
In our calculations, we used a cut-off N max — 8, but we also checked that increasing the number of basis states did 
not change substantially our results for the Mn-Mn separations discussed. The big advantage of using Hermitian 
polynomials as basis states is that for these states one is able to evaluate all matrix elements of the Hamiltonian 
analytically. 

So far we did not discuss the value of the parameters V c , r c , J p d, and r p d, characterizing the the central cell 
correction and exchange interaction. While these phenomenological parameters are not fully determined, they are 
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FIG. 1: Eight lowest lying energy states of the hole localized on the two Mn ions as a function of the distance R between the 
ions. Both Mn spins are pointing in the z-direction. The binding energy of the Mn acceptor (-0.111 eV) is indicated by the 
dashed line. Below R — 10 A the spectrum becomes more complicated then for R > 10 A. 



not completely arbitrary either: they must be chosen in such a way that the model reproduces the experimentally 
observed binding energy of a hole, Emu ~ HI meV^^iiS, and the spin splitting between the bound hole states 
generated by the core Mn spin, -E sp iit ~ 12 meV^S. Therefore, we first performed calculations for a single Mn 
ion, and finally chose a combination of parameters, J p d = 3.7 eV A 3 , r pc i — 0.125 A, V c — 2668 meV and r c = 2.23 A, 
which reproduced the above low energies correctly. 

Having these parameters in hand, we returned to the problem of an Mn2 dimer and computed the spectrum of 
the 'molecule' as a function of the separation R between the two Mn impurities and their spin orientations. Our 
calculations were performed for R in the range 5 Aj R j 20 A. Note that for GaAs lattice constant a is approximately 
a = 5.6 Aand the nearest separation between substitutional atoms in face centered cubic lattice of about 4 A. For 
the sake of simplicity, we performed computations only for the simplest case where the positions of the two Mn ions 
were aligned along the z-direction. In general, however, the spectrum and the effective interaction of the two Mn ions 
depends on the orientation of the Mn bond as well as the orientation of the spins. 

Fig. shows the evolution of the eight lowest-lying states of the Mn2 dimer for the case when both spins point 
in the z-direction. In the absence of the Mn core spin, the ground state of an isolated Mn ion would be fourfold 
degenerat e) 22 ' 25 . This degeneracy is slightly split by a presence of the spin: A classically treated Mn core spin lifts 
this fourfold degeneracy, and results in a small splitting of these four states, which are still well-separated from the 
rest of the spectrum of the ion. 

If we now take approach two Mn ions to each other, then these four states of the two ions hybridize, and give 
rise to eight (bonding and anti-bonding) molecular orbitals. For small enough concentrations it is enough to keep 
only these eight states to build an effective impurity band Hamiltonian for Gai-^Mn^As^. These eight states are 
well-separated from the rest of the spectrum for R > 8 A, but at around R 8 A level crossings occur, and the whole 
envelop function approximation breaks dow for separations smaller than this. 

We can also determine the spectrum of the dimer as a function of the orientation of the two Mn spins, S\ and S2, 
and their separation R. This spin orientation-dependent spectrum is very useful to estimate the exchange energy and 
the anisotropy energy between the two Mn ions. In particular, the spin orientation-dependence of the energy of the 
lowest-energy bonding state just gives the effective interaction energy of the two spins provided that there is a single 
hole on the dimer (roughly corresponding to a hole fraction / = 0.5 %). 

We determined the spectrum of the molecule for two different ('exchange' and 'anisotropy') types of spin configu- 
rations shown in Fig. [21 The corresponding binding energies of a single hole are shown in Fig. [3] For a separation 
R ~ 10 A, the anisotropy energy is very large (in the range of ~ 100 K) but is only about 20 % of the exchange 
energy, which tends to align the two spins ferromagnetically, 8 = 0°. The anisotropy energy prefers an orientation 
perpendicular to the bond, = 90°, in agreement with the RKKY results of Ref. ll6L but in disagreement with the 
results of Fiete et al~, where an 'easy axis' anisotropy has been found. (The origin of this disagreement between the 
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Exchange configuration Anisotropy configuration 



FIG. 2: Configurations used to determine the anisotropy energy and the exchange energy. In the 'exchange configuration' 
one Mn-spin (Si) is pointing in the z-direction and the second Manganese spin (£2) is rotated into the — z-direction. In the 
'anisotropy configuration' both spins Si and S2 are parallel and rotated by the same angle 9. In our calculations spins were 
rotating in the YZ plane. 



results of Ref. |2J and the present computations is unclear). For larger separations the anisotropy energy does not 
decrease too much, but the exchange energy drops by about a factor of 10 for separations ~ 20 A. As a result, the 
anisotropy plays a crucial role for Mn-Mn separations R > 13 A. 
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FIG. 3: Dependence of the binding energy of a hole for two separations (7? = 11 A and R — 19 A) as a function of angle 8. 
Circles stand for anisotropy configuration, squares belong to the 'exchange configuration'. For large separations the anisotropy 
energy is comparable or larger than the exchange energy, while for small separations the exchange energy dominates. Anisotropy 
always prefers an Mn spin orientation perpendicular to the bond, 6 = 90. 

Similar to Refs. IT1IIT2I we can define the exchange and anisotropy energies as 

-Ecxch = ^TT _ E ii ! ^anis = ^TT _ > ( 5 ) 

with the arrows indicating the spin direction with respect to the z-axis. These energies are plotted in Fig. g| Note 
that for separations R > 13 A (corresponding to an active Mn concentration x ac ti V c < 2.2 %) the anisotropy energy 
becomes comparable to the exchange energy. In this range orientational frustration effects discussed in Refs. and 
I23I are expected to be important. 
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FIG. 4: Exchange energy (i? exc h, squares) and non-collinear anisotropy energy (i? an i s , circles) as a function of distance R 
between the two Mn ions. Anisotropy energy is positive and exchange energy is negative indicating that Mn spins tend to line 
up ferromagnetically in the plane perpendicular to the direction joining them. Results obtained using the effective model are 
also plotted as lines. Full lines corresponds to G—4 meV in the effective model and dashed lines to G=10 meV. Anisotropy 
turns out to be the same for both choices of G so the second dashed line is not seen. The effective model give a reasonable 
description of the interaction only for separations R > f6 A. 



III. EFFECTIVE HAMILTONIAN 



As we already mentioned in the previous Section, in the very dilute limit the eight lowest lying states of a pair of 
Mn ions can be constructed from the four lowest-energy acceptor states of the Mn ions, since these are well separated 
from the rest of the spectrum^. Furthermore, the exchange coupling being small compared to the separation between 
the above-mentioned four states and higher-energy excitations, one may attempt to replace it by a local interaction, 
which we write in the second quantized form as 

j pd . ■ ■ -» J2 G Si • = E E G Si • ( c l^c iv ) , (6) 

i=l,2 i=l,2 n,v 

where F{ denotes the effective spin of the bound hole on ion i, and G « 5 meV is an effective coupling that has 
been determined from infrared spectroscopy^. The matrices above are just spin 3/2 matrices spanning the 
four-dimensional Hilbert space of the lowest-lying acceptor states of each Mn ion, and the operators c| create a hole 
on these four states at ion i = with spin component F z = /i. 

Let us set the exchange coupling J p d to zero for a moment. Then, in the spirit of tight binding approximation, the 
most general Hamiltonian can be written within the subspace of these eight acceptor states as 

<-m» = E w^)( c L^+h. c .) (7) 

ju,i/=1...4 
i=l,2 

Note that the hopping terms are not diagonal in F z as a consequence of spin-orbit coupling. Furthermore, spin- 
orbit coupling also generates anisotropy terms, K^u which, however, turn out to be relatively small and can usually 
be neglected. All parameters in Eq. (JJJ depend on the relative position R of the two Mn ions. 

Within the spherical approximation used in Ref. I2H. the angular dependence of and is trivially given by 
spin 3/2 rotation matrices, and the effective Hamiltonian Eq. (J7J) becomes just a function of four parameters that 
depend only on the distance R between the two Mn ions. In the six band model, on the other hand, the numerous 
parameters of the effective Hamiltonian are typically complicated functions of R. However, if the positions of the 
two Mn ions are aligned along the z-direction, than the effective Hamiltonian simplifies a lot due to time reversal 
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FIG. 5: Lowest lying energy states of a manganese pair (lines with symbols) as a function of distance between manganese atoms 
R in the absence of the core spins, e.g., J v d = 0. Dashed lines indicate the results of Ref. 23| Each line is two-fold degenerate. 



symmetry and the C4 symmetry of the Hamiltonian (remember that we set J p d = for the time being), and the 
effective Hamiltonian simplifies as 

H e £-Mn = E t u (R)(c{ tV c 2}U + h.c.) (8) 

u=l..A 

+ YyK{R)^-\)+E{R)]cl vCi , v . 

To determine the four parameters, £1/2, £3/21 K and E, we performed a variational calculation with the previously 
obtained parameters V c , r c , and r p d, but setting J pc i = 0. The results are shown in Fig. EI where for comparison, we 
also show the results obtained within the spherical approximation of Ref. I23I 

While the two approximations give a qualitatively similar spectrum at large separations, the spherical approximation 
badly fails at short separations. There within the six-band model calculation we obtain a bonding and an anti-bonding 
state, each of which has an approximate fourfold spin degeneracy. In other words, at short separations the spin-orbit 
coupling is not very important. In contrast, the spherical approximation gives a large spin-orbit splitting of the states 
even for small separations. 

This difference is easy to understand: for small Mn-Mn separations, the bound states are composed from valence 
band states of all energies. However, while the spherical approximation gives a reasonable approximation for the band 
structure in the close vicinity of the T point, it badly fails at these high energies, where it overestimates the effect of 
spin-orbit coupling. At large separations, on the other hand, the effective Hamiltonian is determined by the tails of the 
wave functions, which are, in turn, composed from small momentum valence band excitations. The structure of these 
states at large distances is therefore roughly captured by spherical approximation, though the results quantitatively 
differ. These observations parallel the ones of Ref. IT2L where we have shown using the RKKY approximation that 
spin-orbit coupling becomes important only for large Mn-Mn separations within that approach too. 

The extracted effective parameters E, K, tn and £32 are shown in Fig. EI The largest parameter is E(R), which 
asymptotically approaches the single ion binding energy. The approach to this asymptotic value is dominated by the 
Coulomb interaction, and E(R) falls off as 

E{R) » £(oc) - ^ . (9) 

The most important effect of spin-orbit coupling is that t\/2 and £3/2 become very different for large separations, 
and \tx/2\ > 1^3/2! implying that holes with spin aligned perpendicular to the bond can hop more easily. This is the 
ultimate reason why spin-orbit coupling favors an Mn spin-orientation perpendicular to the bond. Note that this is 
very different from the result obtained within the spherical approximation, where |ti/2| < 1*3/21) and therefore the Mn 
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spins are preferably aligned along the bond. In fact, the relation | tx /2 1 > K3/2I is also favored by intuition: the weight 
of the valence band orbitals with m = angular momentum in the c£-channel is larger for the F z — ±1/2 states than 
for F z = ±3/2 states. Since the m = orbitals are the ones that point along the Mn-Mn bond, one expects a larger 
hopping in the /i — ±1/2 channels. 
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FIG. 6: Four parameters E (squares), K (circles), t\2 (down triangle) and (up triangle) in two approximations: calculated 
by variational approach with constrains J p d = (lines with symbols) and taken from Ref. [23 (only symbols) . Energy is given 
in Kelvin: 100 K « 8.8 meV. 
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FIG. 7: Energies of the eight lowest lying states of the Mn dimer vs. 8 of a effective model (line) and full calculations (points) 
at R = 11 A. The upper panel corresponds to exchange configuration while the lower panel shows results for the 'anisotropy 
configuration'. 



We shall now proceed and test if the effective Hamiltonian defined by Eqs. © and Q gives indeed a reasonable 
description of the Mn-Mn interaction also in the presence of the Mn spin and for arbitrary spin orientation. To this 
purpose we computed the spectrum of the Mn dimer from the effective Hamiltonian and compared it to the results of 
a full variational calculation with J pc [ ^= 0. The results for separations R — 11 A and R = 19 A are shown in Figs.[7| 
and [SI respectively. Clearly, while for R = 1 1 A the effective Hamiltonian is unable to capture the details of the full 
spectrum, for R — 19 A it gives a very satisfactory description of it. 
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FIG. 8: Energies of the eight lowest lying states of the Mn dimer vs. 9 of a effective model (lines) and full calculations (circles) 
at R = 19 A. The upper panel corresponds to exchange configuration while the lower panel shows results for the 'anisotropy 
configuration'. 



To make a more quantitative comparison, we introduced the average error of the effective Hamiltonian as 

X 2 s (X> - Eff) (10) 

\ — / conf 

where the average is over all spin configurations and the eight lowest-lying states. The quantity \ thus characterizes 
the typical deviation from the real (variationally obtained) spectrum. The obtained value of \ is shown in Table [I] 
The parameters {E 1 K, t 3 / 2 , and t 1 / 2 ) have been extracted from a variational calculation with J pc i = in the way 
discussed above. However, a small error in the largest parameter E results in a large change in x, without changing the 
internal structure and excitation spectrum of the dimer. (The anisotropy or exchange energies, e.g., are completely 
independent of the value of E.) Therefore, to eliminate this error, we optimized the value of E by slightly changing 
it, E — > Sbest- The optimal values of E, and the corresponding x's are shown in Tabled For R > 15 A the effective 
Hamiltonian gives a rather good description, and x 1S only about ~ 10 % of the overall width of the spectrum, 
~ 2|tx/ 2 1 - For R = 11 A, however, the relative error goes above 20 %, and the the effective model's spectrum hardly 
resembles to the one obtained through a full variational calculation with J p d ^ (see Fig.[7J|. 
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-15.9 
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ti/2 [meV] 


-36.9 


-30.4 


-25.6 


-21.6 
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G [meV] 


4.0 


4.0 


4.0 


4.0 


4.0 


E hest [meV] 


-187 


-172 


-159 


-147 


-136 


X [meV] 


15.0 


8.0 


5.8 


4.1 


3.9 



TABLE I: The values of the effective parameters, E, K, ts/2 and ti/ 2 for different R, extracted from a calculation with J p d = 0, 
and the error \ °f the effective model. The typical error \ has been computed after changing the shift E to its optimum value, 
-Ebest, and using G = 4 meV. 

The quality of the effective model approximation can be somewhat further improved by not extracting the values 
of E, K, t 3 /2, and ii/2 from a J p d = variational calculation, but instead, considering them and also G as fitting 
parameters to minimize x- However, even after a full optimization, it is impossible to satisfactorily describe the 
R = 11 A excitation spectrum of the molecule for any spin configuration, while the R > 13 A spectra are not 
substantially improved. 
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IV. CONCLUSION 



In this paper, we studied the 'molecular orbitals' of two substitutional Mn ions as a function of their position and 
spin orientation within the six-band model, using a variational approach. Similar to the results of Ref. Il2t we find 
that spin-orbit coupling effects are important if the separation between the two Mn ions are large (R > 13 A). In 
this regime spin-orbit coupling induces an anisotropy in the spin-spin interaction that prefers to align the Mn spins 
perpendicular to the Mn-Mn dimer's direction. For smaller separations, however, although large, spin-orbit coupling 
induced anisotropy does not seem to be crucial compared to the ferromagnetic exchange interaction. These results 
are in qualitative agreement with those of Refs. 1 1 2ll.il but clearly disagree with those of Ref. 0, where a very small 
anisotropy has b een fou nd. This difference originates from the different cut-off scheme used in Ref. ITU which, as 
discussed in Refs. Illll-j in detail, suppresses back-scattering processes if the cut-off parameter is not chosen carefully 
enough. 

The transition to anisotropy-dominated regions occurs at relatively small active manganese concentrations, x ac ti vo ~ 
2 %. We have to emphasize though that interstitial Mn ions presumably bind to substitutional Mn ions, and therefore 
the active Mn concentration can be substantially reduced compared to the nominal concentration x of Manganes e) 9 ' 10 . 
For an unannealed x — 6 % sample, e.g., with only one third of the Mn ions going to interstitial positions, the active 
Mn concentration can be reduced down to a; ac tive = 2 %, in which regime the anisotropy of the effective spin-spin 
interaction is already very significant. 

Upon annealing, however, interstitial Mn ions presumably diffuse out of the sample, leading to a higher active Mn 
concentration. Therefore, in annealed samples the effective interaction between the Mn core spins should be much 
more isotropic, and the spin of the active Mn ions forms a fully aligned ferromagnetic state. These results are in 
agreement with the experimental data of Ref. |3j, where in some of the samples with small Curie temperature, the 
remanent magnetization could be substantially increased by a relatively small magnetic field, clearly hinting to a 
generically non-collinear magnetic state. We have to mention though that other possible mechanisms have also been 
suggested to cause non-collinear ferromagnetic states^*^. 

We also tested if one can construct a simple effective impurity-band Hamiltonian in terms of spin 3/2 holes hopping 
between the Mn sites to describe Gai_ x Mn x As in the dilute limit, as proposed in Ref. |2j We find that such a 
description only makes sense if Mn-Mn separations are large enough, R > 13 A, and becomes reliable only for active 
Mn concentrations below £ ac tive ~ 1-5 %■ We also determined the parameters of the effective model by performing 
calculations within the framework of the six band model, and found that while earlier calculations done within 
the spherical approximation are qualitatively correct for large Mn separations, they quantitatively fail to reproduce 
the spectrum obtained within the six band model. An interesting result of these six-band model calculations is 
that the effective hopping parameters turn out to be considerably smaller than the ones found within the spherical 
approximation. This implies that the impurity band may survive to somewhat larger concentrations, x ac tivc ~ 
2 — 3 %, as also suggested by angle-resolved photoemission spectroscopy, scanning tunneling microscopy and optical 
conductivity experiments. 

In summary, we find that an effective model description in the spirit of Ref. is only possible if the active Mn 
concentration is less than about a; ac ti V c = 1-5 %. We find furthermore that spin-orbit coupling induced anisotropy 
is extremely large for concentrations active < 2 %, and it likely leads to the appearance of frustrated non-collinear 
states. 



APPENDIX A: SIX-BAND HAMILTONIAN 



For the sake of completeness, we give in this Appendix the effective six-band Hamiltonian used o describe the top of 
the valence band in our computations. Within this approach the valence band holes wave function is a six-component 
spinor, and the so-called Kohn-Luttinger Hamiltonian acts on these six-component spinors^SSiSi 
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In this expression b, c, d, Hh, Hi, and H so denote the following differential operators, 
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and the Kohn-Luttinger Hamiltonian, Eq. (|A1I) . was written in the basis of the following 6 spin states, 
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The first four components of the spinors describe the so-called heavy- and light hole bands on the top of the valence 
band of GaAs. These four bands become degenerate at the V point and have j = 3/2 character. The last two 
components of the spinors describe the spin-orbit split bands. These two bands have a spin j = 1/2 character at the 
T point, where they are separated in GaAs by an amount A so = 0.34 eV from the other four bands. 

The so-called Luttinger parameters, ji, must be determined to give the best agreement with the experimentally 
observed band structure of GaAs. In our calculations we used the set of parameters 71 = 7.65, 72 = 2.41, and 
73 = 3.28, frequently used in the literature^. In the so-called spherical approximation, where only the top four bands 
are kept and the 7j's are set to 71 = 7.65, and 72 = 73 — 2.932^. In this approximation the Hamiltonian acquires an 
SU(2) symmetry. 

For the sake of completeness, let us give here the spin operators of the holes which read in the above basis: 
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